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Abstract 

We investigate the effects of random density fluctuations on neutrino oscillations 
in the Sun environment. We show how the average of certain quantities which can 
be used to describe the MSW effect can be computed analytically. We examine also 
the hypothesis commonly accepted that only perturbations inside the resonance layer 
can have relevance. The average amplitude, which gives the "coherent probability", 
is computed in an analytical and exact way for the case of colored (5-correlated Gaus- 
sian noise: the random perturbation induces a renormalization of the matter density 
which acquires an imaginary part proportional to the fluctuation magnitude in the 
resonance region. Integral equations are given for the density matrix of the system in 
the "optical" approximation. 
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1 Introduction 



In this work we investigate matter-enhanced neutrino flavor transformations, the MSW 
effect ||^, Q , for the case of a exponentially decaying matter density with a colored Gaussian 
noise. 

There are several examples of interesting transport and scattering processes induced 
or modified by the presence of a disordered or random medium: dislocations in crystals 
and the solid-liquid transition; random impurity potentials produce the localization of 
quantum wavefunctions which enables one to understand the transition between insulators 
and conductors; resistance anomalies at low temperature and in presence of magnetic 
fields of "weak localized" electron systems subject to a random potential (|^, ^). The 
same basic localization phenomena explain also in optics the backscattering enhancement 
in electromagnetic wave scattering from a randomly rough surface ([^). Usually, a strong 
dependence on the dimensionality of the problem is observed. In some other cases as in 
nuclear physics the consideration of random matrix hamiltonians allows the simplification 
of otherwise unmanageable systems. 

For neutrino oscillations in presence of random or rapidly twisting magnetic fields 
considerable work has been done already |^, 0, |8|, ^. Random density fluctuations have 
received some attention only recently. 

In Q, a differential equation for the averaged survival probability was derived for 
the case in which the random noise was taken to be a delta-correlated white Gaussian 
distribution. The differential equation was solved numerically and the neutrino evolution 
obtained. Arguments were given which indicate that if the correlation of the matter 
density fluctuations is small compared to the neutrino oscillation length at resonance, one 
obtains the same result as for the case of a delta-correlated noise. In [^] the more realistic 
case of colored noise was considered also in a numerical way and applied to Supernova 
dynamics. An approximate differential equation for the averaged survival probabilities 
was obtained using the hypothesis that the fluctuations should not affect the evolution 
far from the resonance. In |11] , the implications of random perturbations upon the Solar 
neutrino deficit are considered numerically. It is found that the MSW effect is rather stable 
under these fluctuations specially in the small mixing case but in anycase the experimental 
{Am?, cos 9) exclusion curves get modified in an appreciable way. 

In this work we try to develop analytical results for neutrino flavor oscillations induced 
by a random matter density. We are interested mainly in the persistence of the MSW ef- 
fect. The obtention of concrete values for the survival probabilities and phenomenological 
consequences for the solar neutrino problem is left for a subsequent work. The obten- 
tion of these phenomenological consequences is not an easy task , at least, because the 
amplitude of possible density fluctuations in the Sun is poorly known experimentally and 
theoretically. Different arguments can give easily values for it differing by two orders of 
magnitude (|11]): between 0.1% — 10% of the local density. Stronger local inhomogeneities, 
for example at the near-surface dark spots should not be discarded. 

Our basic starting point will be the exact analytical solution obtained in |12] for the 
neutrino oscillation amplitudes in presence of an exponentially decaying matter density. 
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The random component is consider as a perturbation to this solution. Inspired by the 
electromagnetic wave scattering in random media, coherent and incoherent transition 
probabilities are defined. The basic result of this work is that the coherent probability, 
which comes essentially from an averaged amplitude, can be computed exactly in some 
cases. Different integral equations are given for the incoherent probability using an "opti- 
cal" approximation. Approximate solutions valid in a limited range are obtained for them. 
By the very nature of the procedure the results are valid or easi ly generalizable to any 
number of neutrino species. 

The outline of this work is as follows. The first section, after the introduction of some 
known results, is dedicated to explore diverse quantities whose stochastic average can be 
computed exactly or at least by an easy approximation. One of these quantities is the 
determinant of the evolution operator of the system. We use a "naive" argument to show 
the plausibility that the influence of the random perturbation on the neutrino oscillations 
can be described by a complex redefinition of the matter density. This suggestion allows to 
define some ansatz probability. Another quantity is the "total cross section" of the system, 
we suggest the interest in further studying this quantity which can contain information 
on the presence and localization of the MSW resonance. In the same section, the physical 
supposition that only those random perturbations happening in the resonance layer can 
have importance is studied briefly . We argue that this supposition must be taken with 
care, because, as we show, even small phase shifts before the resonance region can have 
some appreciable importance in the flnal survival probability. 

In the next section we define some perturbative expansion for the density matrix of the 
system. The coherent and incoherent parts of this matrix are defined. The coherent part 
being in some sense the zero-order approximation for the full density. Using an "optical" 
approximation, that is discarding a certain class of terms in the perturbative expansion, a 
very general integral equation for the averaged density is derived. In this integral equation, 
basic ingredients are both the coherent density and the averaged amplitude and the two 
point correlation of the matter density perturbation. 

In Section (^) different particular cases are considered. A simpler expression for the pre- 
vious integral equation is given for the case where the random perturbation is 5-correlated. 
A further simpler expression is given for the small mixing case of two neutrino species. In 
this case upper limits for the total averaged probability can be obtained depending on the 
coherent probability and an autocorrelation integral. 

In Section (^) the coherent part is computed. In a first case the small mixing con- 
dition and different approximations are used. Averaging the amplitude amounts to the 
multiplication by a certain slowly time-varying diagonal matrix. In a second particular 
case, it is shown how the average can be computed exactly. The effect of the averaging is 
indeed a complex renormalization of the initial matter density as it was suggested earlier. 
Even if it is a very particular case, it is shown how it can have relevance in more realistic 
computations. Finally, making use of these averages, coherent survival probabilities and 
the "cross-sections" previously defined are computed. 
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2 Preliminary Considerations 
2.1 The non-stochastic solution 

The solution for the neutrino oscihations in solar matter described by the equation 

idtf = {Hq + pq exp{— Xt)uAu~^)i' (1) 

has been given in [p^]. 

For t oo, the solution can be written as (A will be understood generally set to 1) 

i^{t) = ex.p-iHot Uripo) (2) 

For any arbitrary time 

iy{t) = U{t,to)u{0) (3) 

with 

U{t,to) = Us{t)^Us{to), Us{t) = Ur (poe"*) exp-iHot (4) 

The solution can be extended to complex pQ, then Ug is not unitary but still keeps the same 
functional form, the previous expression must be changed (following the same argument 
used to derive Eq.(58) in [^]) to 

U{t,to) = Us{t)-^Us{to) 

In the general case {po real or complex) the coefficients of Ur in the Expression (^) are 
just confluent generalize hypergeometric functions of one order less than the dimension 
of the problem, with argument —ipq/X and parameters which are combinations of the 
eigenvalues of Hq and the elements of the mixing matrix u. For po real, Ur becomes 
unitary and its coefficients adopt a form particularly simple and symmetric. As example, 
we write it explicitly for the two-dimensional case: 

/ 7? Yiz f \ 



Vii 



with the shorthands 



^ exp(z) G* expz F* 

Vu 



G = g{z) = ^iFi(l + Vii /3,2 + /?;z), F = iFi{Vu p,(3;z), 

z = -ipQ, p = Am^/2E 

The algebraic properties of F and G guarantee automatically the unitarity of Ur- Its 
determinant is 

det Ur{po) =exp— ipo (6) 
The survival probability is given by 

Pee(P0,/3,^)= l-2520(l + C20||iFi(i/3C20,l + i/3;-ipo)f) (7) 
or in the small mixing angle limit: 

Pee{po,(3,e)= \\iFi{if3C^e,iP;-ipo)\f (8) 
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2.2 The average determinant 

In this work we are interested to investigate what is the effect of introducing a random 
density perturbation, that means, when p{t) = /3oexp(— At) is changed to p{t){l + 6{t)) 
where 5 is a stochastic Gaussian process of zero mean, characterized completely by a 
certain two point correlation function. 

Having in mind the optical potentials and random matrix approximations used in 
nuclear physics |13] and wave scattering in random media [||], one could think that in our 
case the effect of such introduction should be approximately equivalent to a redefinition of 
the function p{t) which in the most simple case would amount just to a renormalization of 
the constant pQ (to a complex value in general in analogy with the complex wave numbers 
appearing in random media wave scattering). 

We get some hint that such supposition is reasonable if we consider the average value 
of the determinant of the evolution operator of our differential equation. To compute that 
determinant and its average value is a trivial task. 

The evolution operator can be expressed formally as a time-ordered (T) integral: 

U{t,to) =T ex-p-i f dsH{s) (9) 



no 

Its determinant is simply the elementary exponential: 

det U = exp -i [ dstr H{s) (10) 
Jto 

As H{t) is assumed to be a Gaussian process, the statistical average is: 

< det U >=< exp -i J tr H >= exp J < tr-f^ ~\ J J trH{t)trH{s) » 

(11) 

Where «AB >>=< A B > - < A>< B 

So, for the total hamiltonian H defined in Eq.(|^) we have: 

< tr i7(t) >= S + poexp-t (12) 

The constant term S is unimportant, it can be set to zero by a convenient redefinition of 
the zero energy. With 

< 6{t)6{s) >= k{t + s)g{\ t - s \) (13) 

we have 

<< trH{t) trH{s) »= k{t + s)pl exp -{t + s) g{\ t - s |) (14) 

We are interested mainly in the limit t — > oo, we take also to = 0. After some 
elementary integrations we arrive at: 

< det U >= exp — i J dt pQ ^1 — -pQik{t) J g{u)du^ exp(— t) (15) 



^We will use indistinctly the notations < A>, A, or for the average of the variable A 
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We can suppose with generality that g{u) goes to zero quickly for u — > oo , that guar- 
antees that we make a relatively small error if we substitute the upper limit of integration 
by oo. For the case of g{u) being a delta function, this approximation becomes exact. We 
will suppose also that the factor k{t) ~ k, constant. 

So we can interpret that the introduction of the random perturbation induces a renor- 
malization of the initial density which is proportional to the integral of the autocorrelation 
function. 

.kpo 



Po Por = po (^1 - i— giu)duj (16) 

Under this approximation, the average determinant becomes: 

< det U >= exp— zpor (17) 

the renormalization introduces an imaginary part that renders the operator U{por) non- 
unitary, on average; certainly the initial hamiltonian would become non hermitic if we 
would substitute directly it it po by por- Formally we can recover at least the condition 
|< det [/>[=! (but of course not the unitarity condition) adding to the initial hamiltonian 
a diagonal term of the form 



Hr = ipQ exp —t / dug{u) I 
Jo 



I being the identity matrix. This amounts to a (complex) shift in the energy which is 
unobservable. 

We note that for computing Eq.(|^, the unitarity of U has been explicitly used in 



an important way ([12|). Our temptative ansatz is to suppose that the introduction of 
a stochastic term is equivalent to the consideration of a non-stochastic equation with a 
redefined initial maximal density, as obtained before, plus appropriated " counterterms" 
that render the full hamiltonian hermitic. The physical information, the transition proba- 
bilities, can be computed by analytical continuation of the original transition probabilities 
corresponding to the initial equation with the new redefined parameters. 

So under this ansatz the averaged transition probability in presence of the random 
term is (for two dimensions, Eq.(|^)) Pee{por)- 

The rest of this article will be devoted essentially to a more rigorous justification of 
this prescription. We will show that at least for the, so called, coherent probability this 
assumption is true in a particular case. 



2.3 Cross sections and the Optical Theorem. 

As we will see later, while it is relatively easy to compute the average values < Uij >, 
the averaged probabilities, which depend on quadratic quantities <| Uij p>, are nearly 
inaccesibles analytically. 

In fact, apart from survival probabilities, we are also simply interested to study whether 
the MSW effect survives or to which extent gets modified with the introduction of a 
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random perturbation. It would be important to find a quantity that both: gives us some 
information about the existence, amplitude or position in parameter space of the MSW 
effect and its stochastic average is easy to compute . In this sense, we propose to use the 
"scattering" matrix T defined by 

^ oo) = 1 + r (18) 

where Ui is the evolution matrix defined before in some appropriate interaction represen- 
tation. 

T satisfies an " Optical Theorem" : 

UiUj = l = (l + T)(l + r"f) (19) 
T + T^ = -TT'^ (20) 

In particular we define the two quantities: 

a\ = -2mi = E I ^ifc I' (21) 

k 

(72 = -2mrr = 5] I r,fc p=|| T ||2 (22) 

ik 

We can identify the first expression with a total cross section in particle physics (sum 
over all final channels). The second is the sum of all total cross sections over initial 
channels. We expect that the quantities (Ti,(T2 can give us interesting information about 
the scattering process, hopefully the MSW resonance should manifest itself in them. The 
important thing is that, because both depend linearly on T (U), their statistical average 
is rather easy to compute. 

In the mass basis 

= -2^{Uii - 1) = 2(1 - m^i) (23) 

In the weak basis Ui = uUrU~^, and we have a similar expression. The expression for crj 
is particularly simple 

"^al = -2^{{uUrU-^)^i - 1) = 2(1 - mrUrV) (24) 

On the other hand, a2 is basis invariant and (d is the dimension of the problem) 

(72 = 2{d - miUr) (25) 

In Fig.(l^) we plot the quantities '^(7^,(72 for particular parameters together with the 
fe survival probability. We see that both show a prominent peak in the resonance region. 
They reproduce the secondary extrema as well . 

For /?^0or/3— >oo, (72~0 this implies 

2 » mrUr = {1 + cos po)^F + sin po^F (26) 
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For the value of po used in the plot: 1 + cos/9o ~ 2,sm/)o ~ 1/5 (a small but finite value). 
We know also that Pee ^| F p< 1, then - 1 and Pee ~ 1. 

On the other hand, for f3 f3res, (^2 = 4,'^(Tj = 2. miUr = 0,mrC/^y = imply 
~ and we obtain the resonance probability Pee ~ 0. 

Later we will show the same type of plot for the averaged quantities. 



2.4 The influence of the resonance layer 

From a physical point of view, and, as is usually assumed, the only perturbations that can 
have an influence on the final transition probabilities are those which happen inside the 
resonance layer. 

We have supposed Gaussian random perturbations; there is a finite probability that in 
any moment through the neutrino path, a strong local fluctuation makes the density term 
similar to the difference AE provoking in this way a resonant level crossing. Numerical 
studies |10|) seem to show that this possibility in fact doesn't happen easily at least 

for small mixing angles. 

But, we want to show that small changes of the neutrino wave function before arriving 
to the resonance region can have an important effect on the final probability. 

Let's suppose that for any circumstance the flavor components of the neutrino wave 
function acquires a relative phase (p at some ti much later than its creation time. At 
infinite its wave function would be 

= Uli^) S Us{to)i^o (27) 
Where the matrix S is defined by the above equation. 



In Section (5.1) we will see that, although the similarity is not complete, indeed the 
effect of random perturbations can be accounted for by the insertion of a certain matrix 
inside the non-random U. 

In Fig. (^) we show the survival probability as a function of the intermediate time ti (or 
r/ro in the figure) and (p. We see that the phase shifts introduced after the resonance layer 
have no or very little effect. The shifts introduced at the beginning or very clearly before 
the resonance region can have a drastic influence on the final probability. The obvious 
conclusion from this is that if random perturbations affect (even slightly) the phase of the 
wave function long before the resonance region, then they can have an appreciable effect 
on the final survival probability. 



3 Formulation of the Main Approach. 

The density operator of any system defined by a hamiltonian H{t) = Hq + W{t) and 
certain initial conditions, is given, in terms of U, the evolution operator, by: 

p{t) = Uit,to)p{to)U{t,to)^ (28) 
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If the potential vanishes for t — > oo, the asymptotic propagator obeys a free Schroedinger 
equation and can be cast in the form: 

p{t — > oo) = exp —iHot Pas exp iHot (29) 

Pas is a time-independent operator to be determined for any particular problem and initial 
conditions. 

Our objective will be to derive an integral equation for the statistical average of the 
density operator p{t) in the t ^ oo limit for the special type of stochastic potential we 
have discussed in the previous section. 

Using perturbation theory around the known solution of the non-stochastic part of 
W{t), V^{t) =< W{t) >= poexp{-t)V, V^ = Vwe have for U: 

U{t, to) = U\t, to) + ^ i7(") {t, to) (30) 

n 

As we have seen before, for this kind of potential we can factorize the operator in 
"creation" and "anhilation" operators: 

U%to) = U^\t)U^ito) (31) 

so 

U{t,to) = C/f (t)E(i,to)C/°(io) (32) 



oo .J 

S(i, to) = l + Y. dr Vi{Tn) . . . V/(ti) (33) 



{dr = dri X . .. X dTn) 



(34) 



where 

Vj{t) = U^{t){-zV{t))U^\t), Vit) = W{t)-V'{t) (35) 
We define equally the density matrix in the interaction representation 

pi{t) = U^.pUf = E(i, to)pi{to)nt, to)^ (36) 

In the limit t ^ oo, C/°(i) — > exp iHot and we identify 

Pas = Pi{t ^ oo) = E(oo, to) Piito) E(oo, ^o)^ (37) 

The average density operator contain a contribution from neutrinos which have scat- 
tered coherently (or specularly in the electromagnetic wave analogy). This contribution 
is obtained averaging the amplitude U: 

< >=< U> po<U^ > (38) 
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The contribution from the incoherent, or diffuse, component is obtained by subtraction. 



< p^"""'^ >=< p > - < p^"^ > (39) 

An the "randomness" information is included in the operator S, so the Eq.(^) holds 
also for the respective statistical averages: < p >/=< pi >■ 
Taking-^ for the initial density matrix 

Po = u{l,0,...,0)u~^ = V 

the macroscopical z^e^^e transition probability is, 

Pi^it)=tr<pit)>V (40) 

Taking t ^ oo and averaging out oscillatory terms 

P^J = T.iPas,iiVii (~ pas,ii for smaU mixing: Vu 1) (41) 



From Eg. (pg) , the statistical average of pi can be decomposed in a sum of terms of 
the form: 



-pi = E^: 



pq 

P,Q 



< Jvi • • • Vi poi J V} V} > (42) 



It is assumed that the stochastic potential is a Gaussian process with zero mean. This 
implies that the average of terms with an odd number of V's is zero. On the other hand 
the average of the terms with a even number of V's can be decomposed in a sum of two-V 
averages extended over all possible combinations. For example 



12 12 



WW = VV VV+ VVVV + VVVV (43) 
The numbers over the calligraphic V's identifies the pairs which are averaged. 



According to arguments developed in 



13, n] the terms which contain "cross" averages. 



for example the last term in Expression (42), can be discarded to a good approximation 



for pair correlation functions which are significantly non-zero only for a relatively short 
time difference (for example for delta- functions). We will suppose that this is always true 
in our case; in the particular case of the coherent part we will not need this approximation 
at all. 

For convenience we will use from now on a "two time" density matrix, the statistical 
average of pi is given by 



piit,s) = pf'^ + 



^There should not be risk of confusion between p, po as density matrices or matter densities. 
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+ J2 V-'^-VVV-'-Vpoi I V-"--V VV---V 

p,q,r,s-' •' 

/I 2 P 2 "I" X "f" 

V... VVV . . . VVV . . . Vpoi / v...v^v v---v^v v---v^ + 



+••• (44) 

The first summand is the coherent part and its computation amounts to sum the series 
<l+ j Vi+ j Vi j Vi + ... >=< Texp j Vi>=T. (45) 



Note that the average can be performed exactly without neglecting " crossed" pair averages, 
so the last relation Eq.(^5|) can be considered exact. Interchanging averaging and time 
ordering: 



< 



Texp / Vi>=T < exp / >= Texp ^ / / << 1^/(51)^7(52) >> dsids2 (46) 



From this expression an effective potential Vj-^'^ can be defined which summarizes the 
influence of the stochastic potential and such that 

S^// = Texp f Vj^"{t')dt' =^ (47) 
Jto 

Vi^^i') ^ \f« yi{i')yM » ds (48) 

We can derive an alternative approximate expression ignoring crossing terms in Eq. (|4^) . 
The propagator computed in such a case will be called the optical propagator in analogy 
with nuclear physics. It is easy to check that the optical propagator satisfies the equations 

S°P*(i, to) = 1 + /* dr y7*(r)S°P*(r, to) (49) 
Jto 

V°P\r) = f dT2 Vl{T)j:°P\t, T2)V/(t2) (50) 

Jto 

In general Vj^\t,°p^ are not neccesarily the same as the V^^^ , Yf'^^ defined previously. 
Also it is not neccesarily more difficult to compute the former than the latter. In the 
important particular case of a (^-correlated potential both coincide. 

We can compute all the further terms in the expansion of /Tf (Eq.(^4|)) making use of 
the following expression; for an arbitrary non-stochastic operator K, the sum of the series: 

5(t, 5) =< K{t, 5) + r V{t)K(t, 5) + r y(Ti) r V{T2)KiT2, 5) + . . . > (51) 

Jto Jto Jto 
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is equivalent to solve the integral equation 

H 

I 

'to 

and to apply the statistical average. The solution is given by 



S{t, s) = K{t, s) + r dT V{t)S{t, s). (52) 
Jto 



S{t, s) =< T exp /V X K{0, s) + f dr (t exp v] > drK{T, s) 
Jto Jto \ Jto J 



(53) 



which, using the properties of the time ordered exponential and putting -ftr(0, s) = 0, is 
equal to 

S{t,to)= I dT<Texpf V>drK{T,s) (54) 
Jto Jt 



Using this formula iteratively in Eq.(^^ we obtain the expression: 

Jto \Jto / 

+ fY.^Sf{t,r)Vi[r) T S'=^^(r, r2)V/(r2)S^^^(T2, to)po/ f f • 
Jto Jto \JtO 

+••• (55) 

The operators to the right sides of pQi are equal to those in the left sides, changing the 
limit of integration and taking the hermitic conjugate. 

As it can be seen by explicit developing, the series given by Eq.(|55|) is equivalent to 
the following integral equation for the total density pf. 

jAt;s) = T.^ff{tM)poi^'"\sM + f [\^ff{t,n)Vi{n) pi{n,T2) Vj{T2)^'^f^{s,r2) 

Jto Jto 

(56) 

and, at least formally, we have solved the problem in the approximation which neglects 
"crossed" terms. 

4 The Integral Equation for (5-correlated noise. 

Let's suppose that the two point correlation matrix has the form: 

< Vi{t)Vi{s) >= -ik{t)p{t)5{t - s)Vi{t) (57) 
or more generally, for any non-stochastic operator K{t,s), 

< Vi{t) K{t, s) Vi{s) > (= VKV) = k{t)6{t - s)Vi{t) K{t, t) Vi{t) (58) 
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The integral equation for the operator density becomes (Eq.(56)) 



Plit) = E'ffit,to)poi^'"Ht,to)+ / dTkiT)j:'ff{t,T)Vl{T)pj{T)VjiT)^'ff\t,T) (59) 

Jto 



Or recalhng the definition for the representation "I", Eqs.(| 



p{t) = p<=oh(^t) + / dT k{T) U{t,T)V{T)p{T)V\T)U\t,T) 



to 



(60) 



The probabihty P^^ of survival is computed using Eqs.d^). We have the expres- 
sion 

PeJit) = P^fit) + fdr k{r)p\r) trU{t,T)VW)yU\t,T)V (61) 

JO 

Using the property V"^ = V, this expression can be written also in the form 



Pelit) = Pee\t) + fdr k{T)p{Tf tvUnit, t) pR{T)Ul{t, t) 

Jo 



(62) 



Where Xpt = VXV . At this point further approximations are neccesary in order to 
simplify the expressions in eqs.( |6l[]6^ ).. 

4.1 The small mixing approximation. 

In the small mixing, two dimensional case, i.e. taking the limit 

such that Xr = VXV Xupo the integral which appears in Eqs.(|llj6|) b ecomes 

p2(r)fc(r) t f/ii(t,r) (64) 



In this limit, the weak and mass basis coincide and 



pM 

ee 
pcoh 

ee 



Pll 

I Uu P 



So 



Pifit) = P^!Ht) + k{T)p\T)P^°\t,T)P^l{T) 

Jo 



(65) 



Eq.(^) becomes exact for vanishing mixing angles. 



A similar reasoning can be used to derive an equivalent expression for PjJ^ : in this case 
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W is a certain matrix such that = 0, W -V = 0. The final result is 

W = Peti*) + [ dr k{T)p\T)P^;\t, t)P^, {t) (66) 
Defining total probabilities 



pM,coh pM,coh I pM,coh 

ee ' e/j. 



we obtain 



P^{t) = P'^^^it) + f'dr k{T)p\T)P''°''{t,T)P^{T) (67) 

Jo 

This equation is specially useful to set an upper limit on P^ or alternatively on the 
magnitude of the "crossed" terms that we discarded. 
We obtain the strict limit 

P^{t) < P^^^it) + t k{t)p{t) (68) 
Jo 

Setting A;(r) < ke^pPt (/3 < 2 , eventually P = 0) and taking t ^ oo it follows 

pM < pcoh ^ (gg) 

2 p 

In principle P^ = 1, a departure from this value signals a breakdown of the validity of 

the "optical" approximation. 

If we denote by "C" the remaining "crossed" terms not incorporated in the optical 
approximation then 

\c (70) 



J _ _ pcoh _|_ 

So 



1 _ pcoh _ _tpL < c (71) 

An actual estimate can be obtained if the function k{T)p'^{T) is of type exponential, 
so only the values for r ^ of the integrand will contribute significantly. We take the 
approximation 

P^\t,r)Pif{T) P-'^(i,0)P,f (e), 
with e some small number. Under this approximation: 

pM(^) = P^\t) (l + kplP^'fieUt)) (72) 

P^{e) should be typically very small and not very dependent on other parameters; 
Note that for f ^ oo /Lt(t) only varies between Ri 0.5 — 1 irrespective of the exponential 
behavior of k{t). 

We see that the coherent density, or coherent probability, plays an important role not 
only in its own right but also in the computation of the global density or probability. 
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5 Computation of the Coherent probabiUty. 

The effective propagator defined in Eq.(|47|) is 

S^//(t,to) = Texp (^^) j'^dt'k{t')p{t')Vj{t') (73) 



So the effective potential becomes 



yi"{t)=[ '' ]Vi{t) (74) 



In the same case, optical propagator reads: 

rt rrl 



S°f*(t,io) = 1+ f dn r dT2Vj{n)J:°P\n,T2)Vj{T2)^"P\T2,to) 

Jto Jto 

= 1 + ^ fdr k{T)p{T)Vl{T)T."^\T,t,) (75) 

The last relation in Eq.(l75|) is just the integral equation for Y^^f^ which proves that 
both are identical in this case. 
Here: 

V'^f{t) = ^pl^M-m{t)V (76) 

For any function /c(t), a possibility is to compute T,^^^ by perturbation theory, consid- 
ering the addition of a "small" term V^^^ {t) to the hamiltonian H{t) of which we know 
already the exact solution. In Section ^ we will follow this procedure for a particular k{t). 
In fact we will be able to sum the full perturbative series. 

Another possibility, somehow less systematic, will be explored in the next section where 
we will consider the small mixing limit. 

5.1 Small Mixing 

Recalling the definition of Us in Eq.(^) and inserting it in Eq.(7(:) we get 



Vf^t) = exp(-2t)fc(t) ( I/-; I' ) (77) 

In the previous matrix we will further suppose we can neglect the off diagonal terms 
containing oscillating terms of the type exp(i£'i). 

v;/'(«).^exp(-2«)tw(IV'' 1-|«..„|0 ™ 
The time ordered integral is then trivial to perform, taking t ^ oo we obtain 

S^^^(oc,to)=f^^P^-*^^) \ (79) 

V exp(f (</>!- 02))/ 
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where 



p{*o)/po 







p{to)/po 



dxxk{log{x)) I Ur^u{pox) 
dxxk{log{x)) 



(80) 
(81) 



For k{log{x)) ~ kx^ , both quantities 4>i, (j)2 are small, of order 0(1). For an arbitrary 
time t, the expressions will be analogous but changing the inferior limit of integration into 
e ~ 0. So we expect only a very slow dependence on t, this can be seen also from the 
presence of a factor exp(— 2t) in Eq.(|7^). 

In the small mixing limit j Ur^u ^ is the survival probability P^e in the absence 
of random perturbation (Eq.(^)). In a typical interesting case (presence of resonance) 
Pee = Pi ~ 1, for /9oe~*° > Pres and small (Pq ~ 0) for poe~^° < pres- In the absence of 
resonance (formally pres — > oo) Pee ~ 1 for all pc^o- Using this behavior we can estimate 
T.^ff{t oo,to = 0): 



S^^^(P0 > Pr 



^'"{PO < Pr 



exp 



£r^(Pi-Po)-^Po 









exp 



/)2 k 



Pi - Po) + 4^(Po - 1) 



'exp^Pi 






exp^(Pi-l), 



(82) 



The coherent part of the macroscopic probability is (Eqs.(|^|l])) 
For a Tf'^^ = diag{A, B) as before, we have 



■)M,coh 
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(83) 



(84) 



6 The special case k(t)=exp t. 
6.1 The Coherent Probabihty 

It is interesting to consider the case where k{t) = kexpt. Then it is possible to further 
continue the exact analytical expressions. Eq.(|57|) becomes: 

< Vi{t)Vi{s) >= -ikpo6{t - s)Vi{t) (85) 
and the effective potential is 
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The effective potential has the same functional form as the original potential ~ exp —t 
but with the constant po redefined. S"^-^-^ can be computed exactly in a simple way.New 
operators Us,Ur,U are defined by: 

C/f (t)S^^^(t,to)f/i'(*o) = U-\t)Us{to) = U{t,to) (87) 
The new operators are just a renormalized version of the old ones U^, U^. 

U{t,to;po) = & 

In the limit t oo, both C/f , Ug, which depend only on the free hamiltonian Hq, are 
identical, so we get the relation 

S-//(t^oo,to) = Us{to)U^Hto) (89) 

= Ur{por)U}{po) (to = 0) (90) 

The asymptotic density (Eqs.(|29|-^)) is given by 

= UsPoUl (91) 

The coherent probability can be computed using Eq.(^) which amounts essentially to 
redefine po in Eq.(|^). 

The total incoherent part 

P^'*" = 1 - (l Ur,U I' + I Ur,l2 (92) 

is not equal to zero because the matrices Ug , Ur are not unitary now. 



t,to;po 1 



ikpo 



exp iHot Ur 



Po{^ ^)e 



7 Some numerical results 

We see that the choice k{t) ~ fciexpt allows for a fully analytical exact computation of 
the coherent probability. There are no physical grounds for this choice, as there are no 
physical grounds for any other selection, for example k{t) ~ ^2), as long as we don't have a 
very detailed knowledge of the Sun structure. If we suppose that random fluctuations are 
only important if they happen in the resonance layer (but see the warning commentary 
in a previous section) the result of both choices should be approximately equivalent if we 
take ki , k2 such that 

ki exp tres ~ ^2 (93) 
or kiPQ K. k2Pres (94) 

where pres^ires are the position and local density of the resonance layer. 

Outside the resonance region, for t — > oo, ki expt » k2 but in anycase the products 
kiex.ptp^{t),k2p'^{t) ^ 0. 
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For tres > t ^ 0, kiexpt = exp(t — tres)f^2', for a distance {t — t^es) equivalent to 
a quarter of the solar radius that means kiexpt ~ A;2/3, for half of the solar radius 
kiexpt ^ k2/W- The differences between both cases are not so strong as it could be 
thought due to the presence of the extra exponential. 

In Fig.(^) we plot the averaged coherent survival probability P^°^ and the total coher- 
ent probability P'^"^. In the top figure we make use of the Expression ( p^ ) for S^-^-^. In 
the bottom one we plot the raw Formula ([8^). In both figures the probability is strongly 
suppressed for neutrinos created near the origin. In the bottom one, the equivalent fiuc- 
tuation level at the resonance region varies for different creation point r/rg, but, in spite 
of the exponential behavior of k{t), this variation is rather modest. For k = 10~^ as used, 
the equivalent fluctuation at the resonance region r/ro ~ 0.5 — 0.6 is only ~ 7% of the 
local density for a neutrino created at r/rg = 1/10 and « 13% for a neutrino created at 
the center. 

In Fig.(0) we plot instead the averaged coherent probabilities computed using the 
corrected formula por = /Oo(l — ikpres/"^) which guarantees at least uniformity at the 
resonance region for different values of k. 

For small mixing angles, P^°^ shows a moderate damping as k increases. The vari- 
ation for the total coherent probability is much stronger. Beyond the resonance region, 
practically all P™'* « P^f. 

As the mixing angle increases the pattern is different, p™'^ practically doesn't suffer 
alteration but P^°^ is strongly diminished. 

For both angles: for k = 10~^, 10~^ , P'^°^ is practically zero in the inner creation 
regions. Here nearly all the eventual survival probability must come from the diffuse, 
incoherent scattering. 

It is also shown the P*^ defined by Eq.(72). For illustration purposes we choose the 
maximal value P^^{e) = 1. Even in this maximal case P*^ is far from unity, this is particu- 
larly evident for the strongest fluctuations. For large regions the "optical" approximation 
hardly differs from the coherent probability and fails to be an appreciable improvement. 
The lower limit for the "crossed" or "non-optical" corrections , although without much 
value numerically, can be interpreted as indicative that it is precisely around the resonance 
region where the "optical" approximation failure is strongest. 

The average of the cross sections defined by Eq.(24-^) is obtained inserting the average 
value of the matrix U, U^^^ , computed using either Ti^^^ from Eq.(^2|) (Fig.(^), top) or 
Ur{por) (the two bottom figures). In both cases the approximated expression p^es — 
j3cos29 has been used. For the smaller k ~ 10"'* (or 1% fiuctuations) the effect of the 
random perturbation is nearly negligible in both cases. The discrepancy in the large /3 
behavior has the same origin as the difference between Fig.(^) (top) and Fig.(^). For 
k ^ 10"^10~2 the effect can be appreciable. The interpretation of these plots in terms 
of concrete survival probabilities is problematic and subject to further study, but what 
it is clear already from them is that the basic behavior, the very existence, position and 
width of a peak corresponding to the resonance layer remains unaltered. There are no 
discrepancies in this region between the two approaches. 
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In Fig.(|6|), we plot the averaged survival probability for different k using the ansatz 
formula for Pee{por) of Section 2^ and por = po{^ — ipresk/2. There is consistency between 
the behavior expressed here and that of P^"^, P'^"^. If we take at face- value this plot and we 
compare it with Fig.(^) we arrive at the conclusion that for neutrinos created much before 
the resonate region the total survival probability comes essentially from the large avalaible 
quantity of "incoherent" probability. As happens in other physical circumstances, the 
presence of " diffuse" scattering disfavors transitions z^e ~^ and produces an enhancement 
of Pee (state localization). In the same figure, the quick increase of the effect of the random 
perturbations with k is somehow surprising; for ~ 10% fluctuations the effect is already 
quite considerable (" reverse" MS W effect in the case of large mixing angle) . 



8 Brief summary of results and further conclusions. 

The most important result of this work is the derivation of analytical exact expressions 
for the average coherent transition probability for a special case of colored (5-correlated 
Gaussian noise. We have shown that in this case the consequence of the presence of the 
noise is a complex renormalization of the matter density. Other approximative expressions 
for more general cases have been developed as well. It has been suggested an enhancement 
of the survival probability when the incoherent probability becomes dominant. 

It has been proposed for the first time to consider new scattering " cross sections" . It 
has been shown how they are able to describe themselves the MSW effect; clearly further 
work has to be done in this respect. The main importance of these quantities is that their 
statistical average is computable in a simple way. 

The general conclusion is that the MSW effect survives the presence of random per- 
turbations at least for small fiuctuation values. There are indications that their effect can 
become important for large fiuctuations. 
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c2=0.99,z=50 

10 I ' ' ' — ' ' ' — ' — 




Figure 1: The "total cross sections" ^crj (short dash) and a2 (longer dash) for a neutrino 
produced at the Sun {po/X = 50, cos^ = 0.99) as a function of /5 = Am?/2EX (see Eqs.(p4[- 



25)). The continuos line is the survival probability (xlO). 
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arg/2pl 



Figure 2: The survival probability as a function of an arbitrary phase shift introduced in 
different positions along the trajectory of a neutrino created near the Sun core ( cos^ = 
0.99, P = 1, 20) . For (3 = 20 the resonance is situated at r/ro ~ 0.5—0.6. The continuos line 
correspond to r/ro = 1-2: well ahead the resonance region. The longer dashed to r/ro = 
0.6: inside the resonance. The largest variation (both shorter dashed lines) correspond to 
r/ro = 0.5,0.3 at the beginning or clearly before the resonance respectively. 
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Figure 3: The Ve coherent survival probabiUty as function of the neutrino creation point. 
(Top) P^°^ (continues) and P'^°^ (dashed) computed using S^-^-^ given by Eq.(82). Bottom, 
the same probabihties computed using Ur{po{i—ikpQ/2)). {k = 10~^, /3 = 20, cos^ 9 = 0.99 
for both figures.) 
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Figure 4: Coherent probabilities p™" (longer dash), Pg" (continuous line) as a func- 
tion of k and the mixing angle. They are computed using Ur{po{l — ikpresf^))- k = 
10-^,10-^10-2 for upper, middle, lower figures respectively. Right figures: large mixing, 
cos^ 9 = 0.70; left, small mixing cos^ 6 = 0.99. In all the cases (3 = 20. It is also depicted 
P^'^ (Eq.(|7|)) (short denser dash) and Cmin = 1 - P'^°^ - kp'^/2 (short less denser dash). 
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Figure 5: ^it{ (line) and o"2 (dashed line) for a neutrino with the same characteristics 
as in Fig.(|^) in presence of a random perturbation. Two top figures: average amplitude 
computed using H'^-ff from Eq.(82). Lower figures: using Ur{por)- k = 10~^, 10~^ (f« 
1%, 5% fluctuations ) respectively left and right. 



25 



r/rO 



Figure 6: The Vf, survival probability as function of the neutrino creation point computed 
using the ansatz of Section Eq.(0) with por = po{^ — ikpres/'^)- Continuos line: non- 
random probability {k = 0). Dashed lines: k = 10~^, 10~^, 10~^ respectively. For both 
figures /? = 20; cos^ 6 = 0.99 (Top), cos^ 9 = 0.70 (Bottom). 
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